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ABSTRACT 

Finite precision measurement factors the Hilbert space of a quantum system into a tensor 
product -Hcoarse®-fffme- This is mathematically equivalent to the partition into system and 
environment which forms the arena for decoherence, so we describe the consequences of 
the inaccessibility of H& ne as scale decoherence. Considering the experimentally important 
case of a harmonic oscillator potential as well as a periodic piecewise constant potential, we 
show that scale decoherence occurs for inhomogeneous potentials and may explain part of 
the decoherence observed in recent and proposed experiments on mesoscopic superpositions 
of quantum states. 
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Mesoscopic superpositions of quantum states were prepared [1,2] and observed to 
decohere [2] in remarkable experiments reported recently. These experimental successes 
coincided with calls for renewed theoretical investigation of decoherence [3,4] and proposals 
for experimental measurement of the density matrix [1,5] of mesoscopic quantum systems. 
Among the motivations for much of this work is the prospect of quantum computing 
[6]: Not only does the power of quantum algorithms depend on quantum superpositions 
[7,8], but both the ion trap [1] and cavity QED [2] systems used in these experiments 
have been proposed for the physical realization of Boolean quantum computers [9] . Even 
these relatively 'clean' realizations of quantum computation, however, will have runtimes 
limited by decoherence [10]. Current proposals for quantum error correction [11] and fault 
tolerant quantum computation [12] aim to circumvent such limits; detailed theoretical and 
experimental understanding of the decoherence of quantum states may indicate how to 
implement these proposals successfully. 

Careful analysis of the complexity of quantum algorithms considers only finite pre- 
cision specification of transition amplitudes [8], reminding us that only finite precision 
measurements for quantum state observables with continuous spectra are physically real- 
istic. In this paper we consider finite precision measurement of the position of a quantum 
particle evolving according to the Schrodinger equation and demonstrate that in the pres- 
ence of an inhomogeneous potential the measured quantum state generically decoheres— 
independently of any coupling to an environment — a phenomenon which we will call 
scale decoherence. The recognition that quantum measurements have only finite pre- 
cision/limited resolution dates back at least to von Neumann [13]. Coarse graining of 
various sorts, including by finite precision measurement [14], has more recently been used 
in the consistent histories approach to quantum mechanics to obtain decoherent sets of 
histories [15]. The scale decoherence we exhibit here in the reduced density matrix is the 
Schrodinger picture version of this phenomenon and should be apparent in experimental 
measurements of quantum systems having inhomogeneous potentials. These are common 
in the experimental arrangements used for precision measurements of quantum systems 
[3] ; the Paul trap used in [1] , for example, constrains the ions within a harmonic oscillator 
potential. 

Our model for the limited resolution of experimental measurements is the simplest 
possible [13]: the continuous spectrum of the operator of interest is partitioned into con- 
gruent 'bins' [16] corresponding to possible finite precision readings of a measuring device. 
For the one dimensional position measurements with which we are concerned here, we take 
the bins to be intervals of length e > and, in keeping with the computer science conno- 
tations of 'finite precision', write x = ye + z for a position i£R, with y G Z labelling the 
bin and < z < e the position inside the bin. Resolving the mathematical difficulties in 
the usual way [17], we take the set of position eigenstates {\x)} as a basis for the Hilbert 
space H of one particle wave functions. Now notice that writing \x) = \ye + z) = \y) <8> \z) 
defines a tensor product decomposition H = H coarse <g> H^ ne . Since we are taking the 
-fffine factor of the Hilbert space to be inaccessible to measurement, the situation is math- 
ematically equivalent to the decomposition of Hilbert space into system and environment 
factors usual in discussions of decoherence (whether it is conceptualized as damping [18], 



2 



Scale decoherence 



David A. Meyer 



measurement [12,19,20], relaxation [21] or dissipation [22]): small scale position functions 
as an unmeasured 'environment' 'coupled' to the finite precision position 'system'. 

An arbitrary vector \ip) G H need not, of course, factor as a tensor product of vectors 
in if coarse and H& ne . Any quantum state represented as a density operator p G H ® H*, 
however, can be traced over H^ ne to obtain a scale e reduced density operator p := Trfi ne p G 
-^coarse ® H* oarae , describing the quantum state as measureable with finite precision. The 
trace over -£/fi ne is exactly an integral over small scales: for a pure state G H, the 
density function is p(x, x') := ip(x)ifj(x') and the scale e reduced density function is 

P{y,y')-= I ip(ye + z)^(y'e + z)dz (1) 
Jo 

for y, y' G Z. Notice that plane waves are tensor products relative to the decomposition 
of H just as are the position eigenstates. In this case the scale e reduced density operator 
has rank 1: 

p(y,y')= [ e ik (ye+z) e -ik(y'e+z) dz = . jHv-v')*!^ 
Jo 

where 1(0) = e; and thus the scale e entropy S := — Trplogp vanishes as it does for the 
original pure state. 

Now consider a superposition of plane waves: ip(x) = d\e %kxX + a2e tk2X . The scale e 
reduced density function (1) is: 

p(y, y') = a 1 a 1 e ifcl ^- y ') e /(0) + a^e^-^'^/^! - k 2 ) 

+ a 2 a 1 e^ k2y - k ' y '^I(k 2 - h) + a 2 a 2 e tk2(y - y '^I(0), 

where ^ 

I(k 1 -k 2 ) := / e l{kl ~ k2)z dz. (2) 
Jo 

It is clear that in this case, for k± ^ k 2 , p has rank 2. In fact, p is easy to diagonalize: Its 
eigenvectors must be linear combinations <f>(y') = h\e lkiy e + b 2 e lk2y e ; setting 

J2~p(y,y'My') = M(y) 

V' 

we find that for to be an eigenvector A must be a root of the characteristic equation 

= det ( aiai-A aia 2 /(fci - k 2 )/I(0) 

\a 2 dil(k 2 - ki)/I(0) a 2 d 2 - A 

Now the entropy of the scale e reduced density matrix does not vanish: we can solve (3) for 
A and examine the dependence of the eigenvalues (and thence the entropy) on e. From (2) 
it is clear that the relevant quantity is \k± — k 2 \e: when \ki — k 2 \e <^ 1 the scale e entropy 
almost vanishes; only when e is at least comparable to \ki — k 2 \ _1 , which we can consider 
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to be the characteristic spatial scale of ip(x), is there substantial entropy. As e — > oo the 
scale e entropy (in bits) saturates at 1. 

Notice, furthermore, that the eigenvalues A depend on the Fourier coefficients a n 
in the superposition only through their squared norms a n d n . In fact, this is true for 
superpositions of any number of plane waves. For a homogeneous potential plane waves 
are stationary states, so since the time evolution of each is by a phase e~ lU3nt , the scale e 
entropy, even if it is nonzero, must be constant in time; no initial state scale decoheres. 

For an inhomogeneous potential the situation is more interesting: the eigenstates are 
still stationary, but they are not plane waves. Consequently the scale e reduced density 
operator for a wave function expressed as a superposition of eigenstates is more complicated 
to diagonalize. Generically the eigenvalues depend on the phases of the superposition 
coefficients as well as their squared norms and hence the scale e entropy is time dependent. 

Let us examine the case of the harmonic oscillator potential which is present in the 
experimental arrangements described in [1,23,24] and proposed in [5], for example. An 
initial eigenstate does not decohere, of course, since the density operator in this case is time 
independent. So consider an initial superposition of two eigenstates Ui of the harmonic 
oscillator Schrodinger operator, as was prepared in the experiments reported in [23,24]: 
ip(x) = (ui(x) — it3 (a;)) /a/2) say. The scale e reduced density function p(y, y') for such an 
superposition is a complicated and not particularly enlightening expression involving the 
error function, Erf(-). More informative is direct numerical evaluation of the eigenvalues 
and the scale e entropy as a function of time. 

To do so we normalize the potential to V(x) = x 2 /2 and calculate the reduced density 
matrix for e = 1/2. (In order to preserve the symmetry of the problem we shift the domain 
of integration in (1) to [— e/2, e/2).) For numerical purposes the wave function vanishes 
already for \x\ greater than about 5 as is indicated by the probability distributions plotted 
in Figure 1; although e = 1/2 is a significant fraction of the size of the support of the wave 
function it is still small relative to the characteristic spatial scales, approximately 4 and 2, 
respectively, of u\ and U3, so we do not expect the scale e entropy to be saturated. 

Figure 2 shows the scale e entropy (in bits) computed for the time evolution of the 
initial superposition. The scale e entropy initially increases — the state as measured with 
finite precision decoheres. In the longer run, however, the scale e entropy is oscillatory. In 
fact, since the eigenstates in the superposition have periods 47r/3 and 4/t/7, respectively, 
the evolution must be exactly periodic with period tt, as shown in Figure 2. This type 
of oscillatory behavior for the entropy is familiar from studies of the decoherence of a 
finite dimensional quantum system coupled to a finite dimensional environment, e.g., one 
spin coupled to another [20,4]. Here the 'environment' H^ ne is infinite dimensional, but 
only a small finite dimensional subspace — determined by the Fourier components of the 
eigenstates superposed in the initial state — has nonnegligible 'coupling' to the essentially 
finite dimensional 'system' H coarse . This means, in addition, that the scale e entropy is 
bounded above. The maxima of the entropy curve occur at times when the wave function 
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Figure 1. The position probability distributions 
computed from the wave function at time (solid 
curve) and time tt/2 (dotted curve). 



Figure 2. Entropy as a function of time for the 
scale e reduced density matrix of an initial super- 
position of harmonic oscillator eigenstates. 



has the smallest characteristic spatial scale. Figure 1 shows \ip(x, t)\ 2 for t = and t = 7r/2, 
at which times the scale e entropy is minimal and maximal, respectively. 

The minima of the entropy curve in Figure 2 do not give vanishing entropy; that the 
entropy does not vanish at t = (and would be even larger than it is were e > 1/2) is 
uncharacteristic of typical analyses of decoherence where the initial state is often taken to 
factor relative to the tensor product decomposition of the Hilbert space [18-22,3,25]. To 
illustrate this situation for scale decoherence let us consider the case of a periodic potential 
where we can naturally take the initial state to be a plane wave which, as we noted earlier, 
is a tensor product relative to the scale decomposition of H. 

Specifically, we consider the Schrodinger equation defined on the circle S 1 = R/2nZ 
with the potential 



y( x ) = | Vq -tt/2 < x < tt/2 (mod 2tt); 



otherwise. 



(4) 



An eigenstate of energy E is, of course, a linear combination of e l ^ x and e 
and a linear combination of e l ^ E ~ v ° x and e - l ^ E ~ v o x where V(x) 



V(x) 



Ex where 
Vq. The 

usual requirement that the wave function and its first derivative be continuous determines 
the ratios of the coefficients in these superpositions. Requiring the wave function to be 
periodic makes the spectrum of eigenenergies discrete: for Vq = 0, E n = n 2 for n G Z>o; 
for Vq > the eigenvalues split (for n > 0) and are roots of a transcendental equation. 
Let ipn (and i/jq) denote the corresponding symmetric/antisymmetric eigenstates. 

When Vq > a plane wave e lkx for k e Z, although periodic, is not an eigenstate. 
Rather it is a superposition of eigenstates with dominant contributions coming from ones 
with energies near k 2 and near k 2 + Vq. Just as in our earlier example of a superposition 
of two harmonic oscillator eigenstates, this means that the initial plane wave state will 
decohere. To make the scales comparable, consider the case k = 1 in the potential (4) with 
Vq = 3. We expect the dominant components of the superposition to be eigenstates with 
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Figure 3. Entropy as a function of time for the 
scale e reduced density matrix of an initial k = 1 
plane wave in the potential (4) with Vq = 3. 

energies near 1 and near 4. In fact, 
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Figure 4. Entropy as a function of time for the 
scale e reduced density matrix of an initial k = 1 
plane wave in the potential (4) with Vq = 15. 



e lx « - 0.41$j" + 0.55^+ + 0.66^" + 0.18^+ + 0.24^" 

+ O.OlV'a" + 0.03^3" - 0.02^4 , ^ 

where all components with amplitudes of norm at least 0.0001 are included and the en- 
ergies of the first 5 eigenstates listed are approximately 0.52, 3.27, 2.01, 5.41 and 5.84, 
respectively, all of which are reasonably close to 1 and/or 4. The corresponding character- 
istic spatial scales are near 2ir and 7r, respectively, so we set the precision e = n/4 to be 
comparable to the e = 1/2 precision we used in the harmonic oscillator calculation. 

Figure 3 shows the scale e entropy (in bits) for the k = 1 initial plane wave (5). Since 
this initial state does factor as a tensor product relative to the scale decomposition of H, 
the scale e entropy vanishes at t = 0. Then, just as in Figure 2, the scale e entropy increases 
(to a similar maximum value) before decreasing. In this case, however, the multiple higher 
frequency components in the original superposition manifest themselves in the irregularity 
of the entropy curve and prevent the entropy from decreasing all the way to its initial 
(zero) value before it again increases. 

Measuring position with greater precision, i.e., decreasing e, reduces the amount of 
scale e decoherence. On the other hand, increasing the amplitudes of higher frequency 
components amplifies the scale e decoherence. Figure 4 shows the scale e entropy for the 
potential (4) with Vq = 15 but still with e = 7r/4 and initial state a plane wave with k = 1 
as in Figure 3. Now the dominant components are eigenstates with energies near 1 and 
16, and with characteristic spatial scales near 2tt and tt/2, respectively. Including only 
components of the same relevance as in (5) we have now 

e lx w - 0.49^ + O.llV'i" + 0.54zV>r + 0.44^ + 0.29^2 

+ 0.25^ + 0.36^3 - 0.01^4 - 0.01V£ - 0.06^ (6) 
+ 0.02?/>g + 0.02<i/>f - O.OIV^, 
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where the energies of the first 7 eigenstates listed are approximately 0.74, 6.47, 2.92, 15.28, 
11.15, 17.49 and 18.47, respectively, all of which are near the range [1, 16]. Now the scale 
e entropy increases even more rapidly from than in Figure 3, and the presence of addi- 
tional high frequency components makes the entropy curve so irregular as to obscure any 
approximate periodicity in time. The common period of the eigenstates with appreciable 
amplitudes in (6) is so large that the initial state undergoes essentially irreversible scale e 
decoherence. 

These last two computations demonstrate that initial tensor product states scale deco- 
here with the same 'ultraviolet burst' that is familiar from other decoherence calculations 
[3]. Moreover, the last computation shows that when the finite precision 'system' 'cou- 
ples' strongly to multiple dimensions of the small scale 'environment', scale decoherence 
also approaches the irreversible behavior of decoherence by an infinite dimensional thermal 
environment [18-22]. 

Most importantly, we have shown that scale decoherence occurs whenever the poten- 
tial is inhomogeneous, which includes the experimentally relevant situation of a harmonic 
oscillator,* and is mathematically identical to any other form of decoherence. Thus scale 
decoherence is a real experimental effect for finite precision measurements. Although the 
measurements described in [1] are not direct position measurements, the demonstration 
of mesoscopic superpositions depends on the measurement of spin down probability as a 
function of phase separation 4> e [— it, it). At least part of the unexplained decoherence [24] 
in the measurements for initially superposed harmonic oscillator coherent states \a) and 
| — a) with a 2.97(6) may be scale decoherence due to finite precision measurement of 
the relative phase angle 4> — this is consistent with the recognized fluctuations in 4> [1] and 
the finer scale variation of the signal (as a function of (p) in this case than for smaller values 
of a. Similarly, we may expect scale decoherence in the larger a cavity QED experiments 
proposed in [2]. 

Finally, we remark that scale decoherence affects the measurement and phase space 
tomographic reconstruction of the density matrix (or equivalently, the Wigner function) 
for a quantum state [5,23,24]. While not originally conceptualized as decoherence, it is an 
old result [26], revisited in the consistent histories program [27], that smearing the Wigner 
function over a suitable volume of phase space produces a positive probability distribution. 
Finite precision measurement effects exactly such a smearing, although for sufficiently 
small e the resulting scale e decoherence does not eliminate all negative values of the 
Wigner function — and the limited precision need not originate in the quantum mechanical 
uncertainty of simultaneous position and momentum measurements. Since reconstruction 
of the density matrix/ Wigner function is now often performed in quantum optics [28] as 
well as in atomic physics, it may prove useful to consider optical homodyne/phase space 
tomography based on incomplete data [29] from the perspective of scale decoherence. 



Eigenstates, as we have seen, do not scale decohere. Nor, in a harmonic oscillator potential, does a 
coherent gaussian wave packet since evolution preserves its shape; a superposition of such packets 
will, however, scale decohere. 
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